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Abstract. In this paper we study entanglement of the reduced density matrix 
of a bipartite quantum system in a random pure state. 

It transpires that this involves the computation of the smallest eigenvalue 
distribution of the fixed trace Laguerre ensemble of N X N random matri- 
ces. We showed that for finite N the smallest eigenvalue distribution may be 
expressed in terms of Jack polynomials. 

Furthermore, based on the exact results, we found, a limiting distribution, 
when the smallest eigenvalue is suitably scaled with N followed by a large N 
limit. Our results turn out to be the same as the smallest eigenvalue distri- 
bution of the classical Laguerre ensembles without the fixed trace constraint. 
This suggests in a broad sense, the global constraint does not influence local 
correlations, at least, in the large N limit. 

Consequently, we have solved an open problem: The determination of the 
smallest eigenvalue distribution of the reduced density matrix — obtained by 
tracing out the environmental degrees of freedom — for a bipartite quantum 
system of unequal dimensions. 



Key Words: Entanglement, random pure state, Laguerre beta ensemble, ex- 
treme value statistics. 
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1. Introduction 

A bipartite quantum system, is a composite system which can be described as 
a product state of two subsystems; one of which we would like to study and the 
other describing the environment such as noise or heat baths. 

Due to its central role in quantum information and quantum computation, which 
is considered as an indispensable resource, entanglement of a bipartite quantum 
system has recently generated a flurry of activities and has been studied extensively 

m 

There are extensive literature on this topic; here we list those which are of 
immediate relevance to our paper, see for examples [16], [7J, [H], [S], [14], [12] 
and [10]. Statistical properties of such random states are also important in the 
characterization of quantum chaotic systems, see [121 HJ [10] and references therein. 

A bipartite quantum system consists of a system (A) and its environment (B). 
Explicitly, we consider a composite system which is described by a (A^M)-dimensional 
Hilbert space H^ NM ^ = ® ■ Let \ef)f =1 and \ef)™ =1 be the complete 

orthogonal basis for the subsystems A and B, respectively. Without loss of gener- 
ality, we assume that N < M. Any quantum state |$) in the Hilbert space %( NM ) 
can be expressed as the linear combination of \ef) <g> \ef) as follows: 

N M 

where the coefficients Xy € C form a rectangular NxM complex matrix X = [Xijj. 

By a random state we understand that Xij are random variables, so that the 
resulting density matrix of the subsystem A, say, is a random matrix. See (1.2) 
below for a definition. 

The composite state |<I>) is fully unentangled or separable if |$) can be expressed 
as a direct product of two states \$> A ) and |$ B ) drawn from the Hilbert space of A 
and B respectively; that is, 

|*) = 

otherwise, it is referred to as an entangled state. 

We say that |$) is a normalized pure state if and only if associated density matrix 
defined by 

p=|$)($| 

satisfies 

tr[p] = l. 

The reduced density matrix of the subsystem A obtained by tracing out the states 
of subsystem B is found to be, after an easy computation [2l [12] , 

M N 

(1.2) PA := tr B [p] = 5>f \p\ef ) - £ W l3 \e?)(ef\, 

3=1 i,j=l 

where are the entries of N x N square matrix 

W := XX* 

with tr W — 1 implied by the normalization condition that tr [p] = 1. The reduce 
density matrix of subsystem B is similarly defined, ps = tr^p]. 

The fact that tr W = 1 implies the fixed trace ensemble to be introduced below. 
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Because XX* and X^X have the same non-negative eigenvalues, it is not difficult 
to see that the reduced density matrices pa and ps have the same set of non- 
negative eigenvalues 

{Ai}£i 

and satisfy the fixed trace condition 

JV 



5> = i- 



Let vf be the eigenvector of the square matrix W corresponding to the eigenvalue 
Xi. Then the density matrix of the subsystem A can be expressed as 



PA 



N 

E 

i=l 



A similar representation holds for ps- The composite state |$) has the well-known 
Schmidt spectral decomposition [2] 



(1.3) 



1$) = 



N 

E 

i=l 



Ai|0®K 



As we have discussed above, a pure state is random if the coefficients Xy are 
random variables. The simplest and most common choice is to take to be 
independent and identically distributed Gaussian random variables, see |12) for 
detailed exposition. Therefore, the eigenvalues of the reduced density matrix pa 
are distributed according to the joint probability density function (j.p.d.f) of the 
N x N complex Wishart matrix whose trace is constrained to unity, this corresponds 
to the case where the "symmetry" parameter /3 is 2. By the method of random 
matrix theory Q3], it is easily shown that the j.p.d.f of eigenvalues is given by, for 
A^O, j = l,2,..,N, 



N 



(1.4) 



Ps(Xi, A2, • • • , \n) — Cn,mS 



where a = M — N- 



■1-^ > 



-|,/3>0. 



In the situations where (3 = 1, 2, 4 we have 
the real, complex and quaternion fixed trace Wishart random matrix ensembles. 

For general (3, the j.p.d.f. in (|1.4[) can be realized by a tri-diagonal real symmetric 
matrices (with random entries) of the form (see [3]) 

subject to tr[L^] = 1, where Bp is a bi-diagonal random matrix of the form 
/ XPM \ 

X{N-1)P XI3M-I3 

X(N-2)P XPM-2/3 



(1.5) 



B, 



\ 



XI3 XpM-(N-l)P / 



and the random variable Xa has the density 

9I-0/2 
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The normalization constant, which can be computed in closed form, is given by, 

N 



(1-6) C N , M = 



r(fMiV) r(i + § 



nf=o 1 r(f(M- i ))r(i + f(iv-i)) 



See [21] for a derivation of Cn,m- 

The j.p.d.f. given in (jl.4[) . for f3 = 2 can be traced to the work of Lloyd and 
Pagels [H] and of Page [16] . 

When M = N, namely, the two subsystems Hilbert spaces are of equal dimen- 
sion, and (3 = 2, the j.p.d.f. is referred to in [T7] as the ensemble of random density 
matrices with respect to the Hilbert-Schmidt metric in the set of all density matrices 
of size N. The cases where (3 — 1 and (3 — A are also important because these de- 
scribe systems with time-reversal invariance and rotational symmetry respectively 

We assume a > — -| , so that M — N > — 1 . This is chosen so that the normal- 
ization constant (1.6) exists. 

The study of the eigenvalues of the reduced density matrix pa = W is crucial 
for the understanding and utilization of entanglement. In principle, all the infor- 
mation about the spectral properties of the subsystem A, including its degree of 
entanglement, is encoded in the j.p.d.f. given by (|1.4|) . For example, one classical 
measure of entanglement, is the von Neumann entropy defined by 

N 

S = -tr p A In p A = Ajln A. t . 

i=l 

The average entropy (S) — with respect to the j.p.d.f. of (1.4) — was computed by 
Page[16] for large N with M > N. It was found that, 

N 

(S) « In TV , 

V 1 2M' 

which shows that a pure subsystem due to its coupling to the environment is more 
or less completely "randomized" . 

Another important measure of entanglement is the smallest eigenvalue 

A min = min(Ai,A 2 ,...,A A ,), 

and its distribution function. 

This provides, in addition to understanding the nature of entanglement, impor- 
tant information about the degree with which the effective dimension of the Hilbert 
space of the subsystem A can be reduced [T^]. Indeed, the average value (A m i n ) of 
the smallest eigenvalue was studied by Znidaric [50] for the case N = M. Based on 
the exact (A m i n ) for a small values of N, Znidaric conjectured that (A min ) = 1/N 3 
for all N in the complex case {fi — 2). This conjecture was proved in [12] for 
N = M, both for the complex (j3 = 2) and the real (/3 = 1) cases. The problem of 
computing the distribution of A min for unequal dimensions (M > N) was posed in 
|12) . which remains open and it is this problem that we address. 

In this paper we compute the distribution of the smallest eigenvalue A m ; n from 
the joint probability density function of (|1.4[) for general (3. In particular, for the 
most interesting cases where (3 = 1,2,4, we will calculate the distribution of the 
smallest eigenvalue at finite N and in a large N limit to be described later. In the 
situation where (3 — 1, we assume M — N is odd. 
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Let us fix the notations to be used throughout the paper. 
Let Qn,m(x) be the probability that A m i n is > x, that is, 

(1.7) Qn,m(x) ■= Prob [A min > x] = Prob [Ai > x, A 2 > x, . . . , Ajy > x] , 
and the probability density function of the smallest eigenvalues is of course 

(1.8) PnM x ) = -Jj- QnM%)- 

In the situations where /3 = 2, a = (M = N), and f3 = l,a = -\ [M = N), 
closed form expressions of Pn,m{x) and the moments of A m ; n were obtained in 
[T2] . We will extend the above results and solve a problem posed in [12], namely, 
the determination of the probability distribution of A m j n in the case of unequal 
dimensions where M > N. 

Our results are summarized as follows. 

In this paper we extend the known results mentioned above to the cases where 

(1.9) ^ = m>0, meN and (3 > 0. 

We will see that Qn,m{x) and Pn,m(x) are expressed in terms of the Jack symmetric 
polynomials. 

Furthermore, after a scaling 

- JL- 

4iV3' 

followed by the limit N — > oo, we show that, for a fixed m, 
(1.10) 

Q(y) ■= kra^QN^i (4^3) = CX P (~^f) o^i (/3/2) (2m//3;yi, •-•,2/m) \ vl = V2 =...=y m = 

and 
(1.11) 

P(y) := Jirn^ -^Pn.m ) = A m , y m e~^ F^ /2) (2m/ 0+2; y u . . . , y m )\ yi=y 

where 

(1.12) A m ,p = 4™(/3/2)^ 2 + 2m+1 F(1 + /?/2) 



T(l + m)r(l + m + P/2) 



Here 



/2m 
0^1 ( — 52/1, 

is a generalized hypergeometric function. 

Therefore the limiting distributions given by (|1.10[> and (|1.11[> are the same as 
the corresponding results for the Laguerre beta-ensemble, obtained in [S][B] and [3]- 
See also the relevant references in [3]. 

In the next section, we introduce a multiple integral from the j.p.d.f. which 
ultimately determines the smallest eigenvalue distribution for j3 > 0. In section |3j 
by using an alternative method of Mehta 1 1 3J we compute the smallest eigenvalue 
distribution for j3 = 2, and a = M — N a non-negative integer, which we later 
specialize to a = 2. The conclusion can be found in section |4] 
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2. Distribution of the smallest eigenvalue. 

For convenience let us initially take the trace to be t, where t > 0. We define a 
function of x and t as follows: 

r ( N \ N o 

(2.1) I(x,t):= * 5>"* )l[ X ? a II \^-^fd N \ 

J[x,oo) N \ i=1 / i=1 l<j<k<N 

where d^A := dAi • • • cIXn- Clearly, the distribution Qn,m( x ) 1S given by 

I(x,l) 



(2.2) Qn,m{x) 



7(0,1)- 



With the parameter t we may now take a Laplace transform of I(x, t) with respect 
to t. 

Clearly, the numerator of (2.2) is obtained by a Laplace transform of I(x,t) 
with respect to t and followed by setting t — 1 after a Laplace inversion. The 
denominator of (2.2) is 1/Cn,m, with Cn,m as in (1.6). 

To proceed, we first take the Laplace transform I(x, t) with respect to t and 
compare the results to those of Laguerre /3-ensemble, which has been well studied. 
Thus 

(2.3) / I{x,t)e- St dt = / e- s ^^T\ X I a FT \*j ~ d N X. 

J •/[*,«,)* T=l l<j<k<N 

After a linear shift followed a scale change, 

2s/, \ , 



Zt 



(\i~x^ or Xi=x+— Zi, i = l,...,N 



we find, 



(2.4) fli Xt t)e--dt=(iy UN e-^J^a 



where 



(2.5) J{x,-y):= [ f\{x + Zi ye-^ TT |^-^d^. 

•'[0,oo)~ i=1 l<j<fe<iV 

An inverse Laplace transform leaves 

(2.6) I(x, t) = ZT 1 [(£) §MN e -^j{^x, | a )] (t). 
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Hence. 



!N,M 



(x) 



ffs§MN Nx 



J(fx,&a)](l) 



7(0,1) 



c- 



,2s) 



£- 



(1) 



J(0,fa)J(l) 
s f™ J(0,|a) 



(1) 



(2.7) 
(2.8) 



r —MN c 



\MN 



J(0, fa) 



1 ^( 2 «a;, 2 



MAr J(0,fa) 



(1) 



(1 - Nx), 



where we have made used the following properties of the inverse Laplace transform: 

j.a — 1 



(2.9) 

and 
(2.10) 



/r 1 [*-"](«) 



t" 



r(a) 



6>(i), 3?(a) > 0, 



ZT 1 [/( S )e CT ](t)=£- 1 [/(*)] (a + t). 



In (2.9) the Heaviside function 9(t) is 1 for t > and 0(i) is 0, for i < 0. 
We should like to mention that the numerator of 



(2.11) 



j(x, fa) 



is an important but difficult multiple integral, while the denominator can be eval- 
uated as a particular Selberg integral. 

Observe that (|2.11l) can be interpreted as the moment of order (fa) of the 
characteristic polynomial of the equivalent tri-diagonal matrix model of Laguerre 
/3-ensembles [?]. In particular, for f3 = 1, 2, 4, J(x, a/3/2) can be evaluated in closed 
form in terms of Laguerre polynomials, see Chapter 22 of |13) . 

We recall for the reader the definition of the generalized hypergeometric function 
with positive parameter v. This arises in an extensive investigation by Kaneko [9] 
on the multi-variable version of the Amoto's generalization of the Selberg integral. 
See also [BJ, [3]. 

The generalized hypergeometric function of m variables is as follows: 
(2.12) 

i(")...r„ iW 



P x q 



F (l/) (ai, ...,a p ',bt,.. . ,& g ;xi, 



oo r r UU) 

V V [ llK '" [a w R C^Hxi x ) 
2^ 2^ n ,(./) r , ,(„) °* (,xi,...,x m j, 

k=0 |«|=& L°lJt ' ' ' [°«Jk 



where the sum over \n\ = k in (I2.12[) is the sum over all partitions (ki, . . . , /c m ) of 
non-negative integers such that 



(2.13) 
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The generalized factorial function [a]« is defined by 
(2-14) W^:=II - 1 ) 

3=1 V 

with 

(2.15) (a) fe = a(a + l)---(a + fc- 1). 

The to— variable function cj?^ (x\ , . . . , ir m ) is a homogeneous symmetric polynomial 
of degree k and it is proportional to the so-called Jack symmetric polynomials, see 
[9]. A classic reference for Jack symmetric polynomials is [18]. 

In the one variable case, where m = 1, Cr (x) — x k , the generalized hypergeo- 
metric function reduces the one- variable hypergeometric function p F q (x). 

For non-negative integer m , Forrester [6] showed that the quotient of (|2 . 1 1 1) can 
be expressed as a terminating generalized hypergeometric function , 
J(x,m) _ _ 
J(0,m) 



(2.16) -77^ = iF^ /2 \-N;2m/P; Xl ,...,x m )\ Xl=X2= ... Xn 



where m = is a nonnegative integer. Using the above expression, we obtain 
from that 

Qn,m{ x ) 

= T^MN^C- 1 [-|^ /2) (~7V; 2m/ p; x u ..., x m )\ X]= _ 2sx/0 \ (1 - Nx) 

a 00 i r /\/l(^/ 2 ) fc 

= r(fMiv)^- £ r r ; ^/2) Z- 1 [c^Xxn- ■ • , * m )k=...^- 2 ,/, ^iv7l](l-^) 

£Z |fe r(|MJV - fc) [2m//3]^ /2) V P 
(2.17) 

l(/3/2) ^1/3/2) / 1m 



fe=0 |re|=fe 



EE (-1) ^;, a „ T - \ r Il 1^/2) fc! (1 W -**)*""-*- 1 ^-**)- 



Here l m = (1, . . . , 1) where the number of variables is m. We emphasize that (|2.17[) 
is a finite sum. 

Therefore, for finite N, we obtained an expression of the smallest eigenvalue in 
terms of Jack polynomials evaluated at (1, 1, 1) that all the to entries are unity. 

Next, we compute the moments of A m ; n and its limiting distribution after scaling 
with N. 
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Moments o/A m i n . 

From the explicit expression of the distribution in (|2.17j) , one can easily compute 
all the moments of A m i n . For the p-th moment, p > 1, we have that, 

P OO POO 

fip := (Af nin ) = - / x p dQ N , M {x) =P I x p ~ l Q NiM (x)dx 
Jo Jo 

(218) =p g ^ ^ 2 y r(§MN)n P + k) [ _ N ^ dj^) 



to\t^ k K P J T^MN + p)NP+k[2m/fii P/2) 



Limit distribution of X mln . 

We now re-scale x so that x — j^j- Consider a n. If m is fixed, that is, M — N 
does not depend on N, a simple computation produces two useful formulas: 



MN r ah 03/2 



\MN-k-l 



N 



vUMN-k) N V 2 

and 

With (|2TTO|) and (j2~2Ul) . we find, 

Q(2/) = a)Too 0jV ^(4^) 
(2.21) 

= £ £ ^^^1^) e -.f (|) fc = e-^I ^/ 2 ) (2m/A , x , , . , ? , m)|yi 

We see that our smallest eigenvalue distribution with the fixed trace constraint, 
namely, (|2.21[) is the same as that in [6J for the Laguerre /? ensemble without the 
fixed trace constraint. Our result seems to indicate that global constraint does not 
influence local eigenvalue distribution at least after suitable scaling followed by a 
N — ^ oo limit. 

Performing a similar computation, we find, 

(2.22) = A mS y m e^ „ff /2) (2m//3 + 2; Vl , . . . , y m )\ yi=m= ,„ =Vm ^, 

where 

(2.23, A„, a= ^ W 2r^ r(1+ ^^> +OT . 

Note that since P(y) — —Q'(y), (I2.22[) follows from the relation between exact 
expressions of P(y) and Q{y) in [BJ. 

In particular, when m — 0, we find Q(y) = e - ' 9 ^. 

When m = 1, we have, 

Q(y) = 2- 1+ ?r(^)e-^yi-h § _ 1 (^). 
Here I p (x) denotes the modified Bessel function of the first kind pQ. 
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3. another Exact expression at j3 = 2 

In this section, we assume that j3 = 2 and a = AI — N is a nonnegative integer. 
For /3 = 1 or 4, similar but more sophisticated results in Chapter 22 of [13] can be 
used to obtain the corresponding results in this section. 

Our main result is stated below: 



(3.1) Q N , M {x)=T{MN)x 



MN-l£-l 



1 



-.MN 



dot 



k,l=0 



Here (x) is a Laguerre polynomial given by 



N 



(3.2) 



4 p) w = E 

3=0 



N + p\ (~x) j 
N-j 



See the standard reference |19) for Laguerre polynomials. 
Again we begin with Eq. (|2.8[) . but with (3 — 2. We find 



Qn.m(x) =T(MN)£- 1 



1 J(sa;,a) 



„M N 



J(0,a) 



(3.3) 



r(A/A0*V ir 1 J(sx ' a) 



J(0,a) 



S A1N 



(1 - 2Vz) 

(1 - JVx). 



J(0,0) 



The integral (Eq.(17.6.5), Page 321, [13]) implies 



N 



(3.4) 

and hence 
(3.5) 



J(0,a) = ]jT(l+j)T(a + j), 

3=1 



N 



.7(iu)i _jj n./) 



iV 

no-)c 



J(0,a) ±±r(a + j) 
On the other hand, it follows from Eq.(22.2.28), Page 416, 13 , that 

d l 



(3.6) 



)aJV J(-a;,a) _ /yj A 



det 



fix 



l C N+k (x) 



a-l 
k,l=0 



where Cj(x) is the monic polynomial related to the Laguerre polynomials as follows: 
(3.7) C ] (x) = (-iy 3 \Lf\x). 

By an elementary differential-difference relation satisfied by the Laguerre poly- 
nomials, 



(3.8) 



-lJf(x) = -I^>(x) 1 



(p+i)/ 



dx 
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we obtain 

. aN J(-x,a) 



J(0,0) 



II A l[(-l) N+l (N + l)ldet\^Lf +k ( X ) 



1=0 



{ _ 1)aN -q 



Q-l 



det 



-(0 



(3.9) 

Note that 
(3.10) 



(-l) aN Y[(l + l) N det 



2=0 



L 



N+k-l( x ) 
(0 



L 



N+k-l 



dx l 

a-l 

k,l=0 

a-l 
k,l=0 



a-l 

k,l=0 



n (■?■)«= nV+ik- 

j=l (=0 

Combining Eqs. (j3~5|) . (|33|) and (I3TTU)) . we obtain, 



(3.11) 



J(—x, a) 



= det 



a-l 



J(0,a) L J,T "- " 'Jfc,i=o 

Furthermore, from Q3.3p . we arrive at the expression, 



Qn,m(z) = r(MiV)£" 



1 



det 



k,l=0- 



(1 - Ax) 



MN . — 1 



T(MN)x MN C 



1 



(3.12) 



T(MiV)a; 



MN — 1 



{sx) MN 
1 



det 
det 



(0 

N+k-l 



(sx) 



a-l 



-(I) 

'N+k-l 



(s) 



k,l=0 
a-l 



k,l=0 



(1 - Mc) 

1 - Nx^ 

x 



Fr 1 ). 



where we have used a property of inverse Laplace transform: 

1 

.6- 

Now we focus on the case a = M — N — 2. Note that since, 

N 



(3.13) 



we have 



4 P) (-*) 



-f(^)](t) = C- 1 [f( s )}(bt). 
- M — N = 2. Note that s 



An 



E 

3=0 



(p + i)3 i! 



- , (zfl! y L^Qj Lf)! N y (-^- !)i (-«)' y 1 (-^+1)3 (-^ 



'• ft ^ 



EE 



i=0 v '* 3=0 v /J t=0 

(-«)<+' (AT + 1)(—N)i(—N)j - N(-N - l)i(-N + I), 



3=0 



(2)3 



(l)*(2)i 



(3.14) 



EE 



(- s )*+3> A r) i (-A0 J -(Ar+l)(l+j-i) 



(1)<(2) 3 - 



» j 

Here we have used the fact 



N+l-i 



(N + i)(-N)i(-N) s - N(—N - l),(-iV + 1), = (-A04-A03 



(N +!){!+ j-i) 
N + l-i 
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In addition, since we are only interested the case of very large N, we can always 
write the sum in the form of (3.14) for fixed i,j. Thus, it follows from Eq. (|3.12j) 
that 



Qn,m{x) 
= T{MN)x MN - 



'EE 



(_1)<+J (-NU-N), (N + 1)(1 + j - 



(l)i(2)j 



N + 1 



1 



(-1)^ (-NU-N) J (N + j- i) 



s MN-i-j 
1 1—Nx \ 



1 - Nx\ 
x J 



(3.15) 

= (1 - Nx) MN - 



X E 



(l)i(2)j W + l-i r(MJV-i-j) 

(-^(-iV^ (TV + 1)(1 + j - i) T{MN) 



1 - Nx\MN-\-i-j 



(l)i(2) 



iV + l-i 



T(MN — i — j) ^_ N ^ 



i+j ' 



Setting x = j^-. Consider a fixed We now deal with the limit of the 



individual factors on the right-hand side of Eq. (|3.15|) as N — > oo. 

y 

i \ ; 



It is easily seen that (with x — t^j), as N — > oo, 



sMJV-l 



(l-Nxy 

(1 - Nx) l+J 
(-NU-N), 

T(MN) 
T(MN -i-j) N 2i + 2 i 



p-V/4 

1, 

(-1)*+', 
1 



Hence, 



lim Qn,m(-ttp>) 

N^too AN 3 



= e ~ v/A E 



(y/4) i+ Jl+j-z 
(1M2), 



-y/4 



E 



(2//4) 



1 



(» - + l)! 



i,j = 

' ^ ^ j!?! 4 e ^ (i - 1)H! ^ 

(3-16) =e^/ 4 (/ 2 (Vy)-/ 1 2 (Vy)), 

where the modified Bessel function of the first kind is given by 

(x/2) 2k 



(y/4)* 



(3.17) 



I p (x) = (x/2yJ2 



k=0 

4. Conclusion 



k\T(p + k + 1) ' 



In this paper we have studied the exact probability distribution of the smallest 
eigenvalue of the density matrix of an entangled random pure state. This turns out 
to be the same as the smallest eigenvalue distribution of the fixed trace Laguerre 
random matrix ensemble. We obtained, for a bipartite quantum system composed 
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of two subsystems whose Hilbcrt spaces have dimensions M and N respectively, the 
exact smallest eigenvalue distribution for all M > N (in the (3 = 1 case we assume 
M — N is an odd integer). The distributions are expressed as the evaluation of a 
certain hypergeometric functions at special points and solve an open problem in 
quantum entanglement. 

Our result not only provides important information on entanglement of a bipar- 
tite quantum system in a random pure state, but also demonstrates the intimate 
relations between entanglement of bipartite quantum systems and the fixed trace 
Laguerre ensemble. We may conclude based on our results that in a broad sense 
the global constraint does not influence local correlations at least in a certain large 
N limit. An indication that this may have wider validity can be found in [10) where 
the authors studied the kernel and the n— point correlation functions of the fixed 
and bounded trace Laguerre ensembles. There it was found that the suitably scaled 
kernel in the large N limit converges to the unrestricted kernel in the bulk, hard 
and soft edges, using the terminology of Tracy and Widom. However, we should 
like to mention that the "matching up" of the kernels in the constrained and the 
unrestricted ensembles does not imply that the "matching up" of the corresponding 
distribution functions (in this instance the smallest eigenvalue distributions) . Only 
after a rather elaborate computation shows that this is the case. 

Although there is currently no obviously effective characterization of the degree 
of entanglement, the von Neumann entropy, however, is considered to be useful 
as a measurement of entanglement in bipartite quantum systems. The distribu- 
tion function of the von Neumann entropy may be successfully tackled using our 
approach. 
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